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The coherent contribution of all neutrons in neutrino nucleus scattering due to the neutral current 
is examined considering the boron solar neutrinos. These neutrinos could potentially become a source 
of background in the future dark matter searches aiming at nucleon cross sections in the region well 
below the 10^^"pb, i.e a few events per ton per year. 
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INTRODUCTION. 

The universe is observed to contain large amounts of dark matter [J, [3| , and its contribution to 
the total energy density is estimated to be ~ 25%. This non-baryonic dark matter component, 
responsible for the growth of cosmological perturbations through gravitational instability, has still 
not been detected directly. Even though there exists firm indirect evidence from the halos of dark 
matter in galaxies and clusters of galaxies it is essential to detect matter directly. 

The possibility of direct detection, however, depends on the nature of the dark matter constituents, 
i.e the WIMPs (weakly interacting massive particles). Supersymmetry naturally provides candidates 
for these constituents [1, 0) @, B 0] ■ most favored scenario of supersymmetry, the lightest 

supersymmetric particle (LSP) can be described as a Majorana fermion, a linear combination of 
the neutral components of the gauginos and higgsinos. Other possibilities also exist, see, e.g. some 
models in universal theories with extra dimensions Q. Since the WIMPs are expected to be very 
massive {mwiMP > 30 GeV) and extremely non-relativistic with average kinetic energy (T) ~ 50 
keV (rniv/ Afp/100 GeV), a WIMP interaction with a nucleus in an underground detector is not 
likely to produce excitation. As a result, WIMPs can be directly detected mainly via the recoil of 
a nucleus {A, Z) in elastic scattering. The event rate for such a process can be computed from the 
following ingredients: 

1. An effective Lagrangian at the elementary particle (quark) level obtained in the framework 
of theprevailing particle theory. For supersymmetry this is achieved as described, e. g., in 
refs. [sl, 0|. For Kaluza Klein theories in universal extra dimensions see, e. g., some recent 
calculations [l3|- Invariably this ingredient is the most important element, but at present, 
unfortunately, with the greatest uncertainty in getting the event rate, especially since the 
WIMP mass is quite uncertain. 

2. A well defined procedure for transforming the amplitude obtained using the previous effective 
Lagrangian from the quark to the nucleon level, i.e. a quark model for the nucleon. This step 
in SUSY models is non-trivial, since the obtained results depend crucially on the content of 
the nucleon in quarks other than u and d. 

3. Knowledge of the relevant nuclear matrix elements [lll.[l^. obtained with reliable many-body 
nuclear wave functions. Fortunately, in the case of the scalar coupling, which is viewed as the 
most important, the situation is a bit simpler, as only the nuclear form factor is needed. 

4. Knowledge of the WIMP density in our vicinity and its velocity distribution. Since the essential 
input here comes from the rotation curves, dark matter candidates other than the LSP are 



also characterized by similar parameters. One does not know for sure what these parameters 
are. The most common models are isothermal Maxwell-Boltzmann (M-B) distributions llSlJl^ 



151 . Il6l | , anisotropic velocity distributions described by Tsallis type functions [17|, [18 




i.e. functions which in some limit lead to the M-B distributions, and variants of the M-B 



distributions arising when dark matter is coupled to dark energy [2l|. Non isothermal models, 
like those arising in the Eddington approach have also been considered [HSHKil. These 
various models give similar (time averaged) rates and differ only in the predictions regarding 
the modulation. So one may assume that the uncertainties in this case are quite small. 

The particle physics in conjunction with the structure of the nucleon provide the nucleon cross 
sections. Since, as we have already mentioned, the particle physics parameters most likely will 
result in very small cross sections, the most ambitious future dark matter experiments like the 
XENON-ZEPLIN aim at detecting 10 events per ton per year. At this level one may encounter 
very bothersome backgrounds. One such background may come from the high energy boron solar 
neutrinos (the other neutrinos are characterized either by too small energy or much lower fluxes). 

During the last years various detectors aiming at detecting recoiling nuclei have been developed in 
connection with dark matter searches [2^ with thresholds in the few keV region. Recently, however, 
it has become feasible to detect neutrinos by measuring the recoiling nucleus and employing gaseous 
detectors with much lower threshold energies [l^. Thus one is able to explore the advantages 
offered by the neutral current interaction, exploring ideas put forward more than a decade ago [28|. 
Furthermore this interaction, through its vector component, can lead to coherence, i.e. an additive 
contribution of all neutrons in the nucleus (the vector contribution of the protons is tiny, so the 
coherence is mainly due to the neutrons of the nucleus). 

In this paper we will derive the differential neutrino nucleus cross section and the associated 
event rate for the elastic (coherent) neutrino-nucleus scattering. Then we will utilize the available 
information regarding the energy spectrum of solar boron neutrinos and estimate the expected 
number of events for light as well as heavy nuclear target. Finally we will compare the recoil 
spectrum and total event rate associated with WIMPs with that due to neutrinos. 



A BRIEF DISCUSSION OF THE RATES FOR DIRECT WIMP DETECTION 

Before proceeding with the evaluation of the event rate for nuclear recoils originating from the 
neutrino nuclear scattering we will briefly discuss the WIMP-nucleus recoiling events. We begin by 
saying that the shape of the differential event rate for WIMP detection cannot be precisely estimated, 
since, as we mentioned in the introduction, the WIMP-nucleon cross section is not known. Especially 
its dependence on the WIMP mass is not known. The non-directional differential rate folded with 
the WIMP velocity distribution is given by [H, 0, HH, [H] : 



(f) - '-^i^VW) f JfL/(,,,,)^,3, (1) 

\ du I ArriN J \/ (v^) du 

where f{v,VE) essentially is the WIMP velocity distribution in the laboratory frame with Ve the 
velocity of the Earth. 

The differential cross section is given by: 

d<j(u, v) == ^^^^ [i^sF{uf + f:,p,„Fn{u)] (2) 

where u the energy transfer Q in dimensionless units given by 

u=-^ , Qo = [rUpAb]-^ = MeV (3) 

Qo 



with b is the nuclear (harmonic oscillator) size parameter. F(u) is the nuclear form factor and Fii{u) 
is the spin response function associated with the isovector channel. 



The scalar cross section is given by: 
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(since the heavy quarks dominate the isovector contribution is negligible), afr is the LSP-nucleon 



scalar cross section. The spin cross section is given by 

-'spin \ / \ 
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Cspin the nuclear spin ME. 

Integrating over the energy transfer u we obtain the event rate for the coherent WIMP-nucleus 
elastic scattering, which is given by [H, Him, [12: 
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In this work we will not be concerned with the spin cross section. The parameter tcoh takes into 
account the folding of the nuclear form factor with the WIMP velocity distribution, hcoh deals with 
the modulation due to the motion of the Earth and a is the phase of the Earth (zero around June 
2nd). 

The number of events in time t due to the scalar interaction, which leads to coherence [2jy], is: 
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Assuming a constant nucleon cross section we get the differential rate indicated by Figs [T] and [2l 
This shape was obtained in the coherent mode, but we expect to have a similar shape due to the 
spin 2^, 33 1 . The essential difference is that one needs as input a much larger nucleon cross section 
due to the spin. Indeed for a heavy target the spin mode has no chance, if the nucleon cross section 
due to the spin is of the same magnitude with that associated with the coherent mode. 
The total (time averaged) coherent event rate is shown in Figs [5] and [H 



ELASTIC NEUTRINO NUCLEON SCATTERING 



The cross section for elastic neutrino nucleon scattering has extensively been studied. It has been 
shown that at low energies it can be simplified and be cast in the form: [33],[3^: 
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FIG. 1: The differential event rate for the coherent process, in arbitrary units, as a function of the parameter 
ay'M, where a — 1.410^''' (/ir(A)c'^b/ {Tic) with /ir(^) the WIMP-nuclcus reduced mass and h the size parameter 
for the nucleus, u is essentially the energy transfer Q. u = Q/Qo, Qo = 4.110^A~*/^. Due to the nuclear 
form factor not all the range of u is exploitable in direct WIMP detection. For '^^'^Xe, e.g., effectively there 
is an upper cut off value indicated by dotted line, fine line and thick line for a WIMP mass of 30, 100 and 
200 GeV respectively. Any lower cut off is due to the threshold. 




Q — > keV 

FIG. 2: We show the differential event rate dRcoh/dQ for the coherent process , as a function of the energy 
transfer, for a WIMP mass, rrix, of 100 GeV in the case of "iXe. 



where niN is the nucleon mass and gv, qa arc the weak coupling constants. Neglecting their 
dependence on the momentum transfer to the nucleon they take the form; 



gv = -2sm^ Ow + 1/2 « 0.04 , 5^ = — , {u,p) 



(11) 
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gv^-l/2, gA^-— , {iy,n) (12) 

In the above expressions for the axial current the renormalization in going from the quark to the 
nucleon level was taken into account. For antineutrinos qa —gA- To set the scale we write: 

^ 5.14 X 10-41 ^ (13) 
2tt MeV ^ ' 



The nucleon energy depends on the neutrino energy and the scattering angle, the angle between the 
direction of the recoiling particle and that of the incident neutrino. In the laboratory frame it is 
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FIG. 6: The same as in Fig. [5] for a detector energy threshold of 10 keV. 



given by: 

J-N = -, —ETT? 7E av? ' < 6^ < 7r/2 (14) 

[mN + Ei,Y - [E^ cos Oy 

(forward scattering). For sufRciently small neutrino energies, the last equation can be simplified as 
follows: 

2(^.cosg)^ 
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The above formula can be generalized to any target and can be written in dimensionless form as 
follows: 

2C0S2 6I _ Trecoil _ 
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In the present calculation we will treat Xi, and y as dynamical variables, in line with CDM recoils. 
One, of course, equally well could have chosen x^, and 9 as relevant variables. 
The maximum energy occurs when 9 = 0, i.e.: 

2 



Umax ~ /I I -, / \2 1 ' {^^) 

(1 + i/xi,y - 1 

in agreement with Eq. (2.5) of ref. [33|. One can invert Eq. [T5l and get the neutrino energy 
associated with a given recoil energy and scattering angle. One finds 
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The minimum neutrino energy for a given recoil energy is given by: 
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in agreement with Eq. (4.2) of ref. [3J|. The last equation is useful in obtaining the differential 
cross section (with respect to the recoil energy) after folding with the neutrino spectrum 



COHERENT NEUTRINO NUCLEUS SCATTERING 

From the above expressions we see that the vector current contribution, which rt iay lead to co- 
herence, is negligible in the case of the protons. Thus the coherent contribution [36| may come 
from the neutrons and is expected to be proportional to the square of the neutron number. The 
neutrino-nucleus scattering can be derived in analogous fashion. It can also be obtained from the 
amplitude of the neutrino nucleon scattering by employing the appropriate kinematics, i.e. those 
involving the elastically scattered nucleus and the substitution 



P P ^ / 2 ^ P' 



with q the nucleon momentum and p the nuclear momentum. Under the above assumptions the 
neutrino-nucleus cross section takes the form: 

+ {My - M^ni - + (Ml - M^)^!|p] (19) 

where My and Ma are the nuclear matrix elements associated with the vector and the axial currents 
respectively and Ta is the energy of the recoiling nucleus. The axial current contribution vanishes 
for 0+ 0+ transitions. Anyway it is negligible in front of the coherent scattering due to neutrons. 
Thus the previous formula is reduced to: 



f da \ GpArriN 
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{N^/A)F,oh{TA,E,), (20) 



with 



F..u{Ta. E.) = F\e) (l + (1 - - (21) 

where F{q^) — F{T\ + 2AmMTA) is the nuclear form factor. 

The nuclear form factor makes a sizable contribution in the case of the more energetic supernova 
neutrinos. In the case of solar neutrinos, due to the small nuclear recoil energy, the form factor is 
expected to play a minor role numerically. 

By setting the form factor equal to unity in Fcoh we obtain fcoh- The latter is shown in Fig. [51 
Note that the maximum recoil energy for boron neutrinos cannot exceed the 4 keV in the case of Xe 
and 17 keV in the case of S. The differential cross section is shown in Figs fT2l and fT3l Note the rapid 
fall of the cross section with recoil energy. The cross section has been almost completely depleted 
beyond energies 1 and 5 keV for a intermediate (^"^^Xe ) and hght ("^^S) targets respectively. 



THE EVENT RATE 



To proceed further we must convolute the cross section with the neutrino spectrum. From the 
neutrinos emitted by the sun only the boron neutrinos have high enough flux with sufficiently 
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FIG. 9: The function fcohiTA), obtained from F^ohiTA) after setting the form factor and the quenching 
factor equal to unity, for A=131 on the left and A=32 on the right as a function of the recoil energy ( for 
neutrino energies 3, 6, 9, 12 and 15 keV increasing to the right). According to Eq. p6p the maximum recoil 
energy is increasing as the neutrino energy increases. 




FIG. 10: The boron solar neutrino spectrum. 




FIG. 11: The boron solar neutrino flux in units of cm s ^ in a Log-Log plot. 



high energy to lead to nuclear recoils, which could become relevant in dark matter searches. The 
normalized boron neutrino spectrum is shown in Fig. 1101 The corresponding flux is shown in Fig. 
[TTl The obtained differential cross section is shown in Figs [12] and [131 

Integrating the differential cross section down to zero threshold we find the event rates given 
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FIG. 12: The neutrino induced differential cross section in units cm^ keV~^ as a function of the recoil energy 
in keV in the case of the target ^"^^Xe. Note that essentially all the contribution comes from the recoil energy 
region < Ta < 1 keV. The effect of the form factor is invisible in the figure. 
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FIG. 13: The same as in Fig. [T2]for the target ^^S. In this case the recoil energy region is wider < Ta < 5 
keV, but still small. 



in tabic |T1 The event rates are almost two orders of magnitude smaller than the rates for WIMP 
detection obtained with a nucleon cross section of lO^^pb. Thus such neutrinos cannot be a serious 
background for WIMP searches in the region 10^^ — lO^^^pb. In any event, as we will see below, 
the neutrino induced recoils are less of a background problem in the realistic case of non zero energy 
threshold. 

Sometimes for experimental purposes one may have to focus on a resticted region of the recoill energy 
spectrum. To be specific let us consider a typical low recoil energy region, e.g. 2 keV< T4 <4 keV. 
Clearly from Fig. 1121 one can see that the neutrino background is very small in this energy region. 



TABLE I: Comparison of the event rates for boron solar neutrino detection with those of WIMP detection 
rates. In evaluating the latter we assumed a nucleon cross section independent of the mass. The kinematics 
were obtained assuming two WIMP masses, namely 100 and 300 GeV. NoFF means that the nuclear form 
factor was neglected. 
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We find that the WIMP event rate in this restricted energy region is substantially reduced, but it is 
not suppressed as much as the neutrino rate (see table |T| . 

We should mention that the obtained rates are independent of the neutrino oscillation parameters, 
since the neutral current events, which we considered in the present calculation, are not affected by 
such oscillations. 



QUENCHING FACTORS AND ENERGY THRESHOLDS 



. The above results refer to an ideal detector operating down to zero energy threshold. For a real 
detector, however, as we have already mentioned, the nuclear recoil events are quenched, especially 
at low energies. The quenching factor for a given detector is the ratio of the signal height for a recoil 
track divided by that of an electron signal height with the same energy. We should not forget that 
the signal heights depend on the velocity and how the signals are extracted experimentally. The 
actual quenching factors must be determined experimentally for each target. In the case of Nal the 
quenching factor is 0.05, while for Ge and Si it is 0.2-0.3. For our purposes it is adequate, to multiply 
the energy scale by an recoil energy dependent quenching factor, Qfac{TA) adequately described by 
the Lidhard theory fST']. More specifically in our estimate of Qu{Ta) we assumed a quenching factor 
of the following empirical form ^SSj-iS?]: 



Qfac{TA) = ri 
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n ~ 0.256 , r2 ~ 0.153 



(22) 



The quenching factors very much depend on the detector type. The quenching factor, exhibited 
in Figs [ni and [T5l for recoil energies of ^^^Xe and ^^^S respectively, were obtained assuming the same 
quenching of the form of Eq. (|22p . In the presence of the quenching factor as given by Eq. ( the 
measured recoil energy is typically reduced by factors of about 3, when compared with the electron 
energy. In other words a threshold energy of electrons of 1 keV becomes 3 keV for nuclear recoils. 
Accordingly, the event rates for neutrino recoils are reduced much, as seen from Figs [16] and [17] 
below. On the other hand the WIMP recoil events are not reduced much, since the recoil energy is 
well above threshold. 

The above rates were obtained in the case of zero threshold. Due to the relatively low recoil energies. 




however, the effect of threshold is crucial (see Figs [TBI and fT7|). One clearly sees that the observed 
events are an order of magnitude down if the energy threshold is 1 keV (2 keV) for Xe(S) respectively. 
Thus with quenching most signals get below the threshold energy of « 1 keV. On the other hand the 
WIMP event rates are almost unaffected, unless the threshold energy becomes larger than 5 keV. 



CONCLUDING REMARKS 



In the present study we considered the elastic scattering of WIMP-nucleus interaction and the 
corresponding elastic scattering of boron solar neutrinos. The former are favored by an enhance- 
ment due to coherence of all nucleons, while the latter by due to the neutron coherence resulting 
from the neutral current interaction. The latter may become a source of background, if the WIMP 
nucleon interaction turns out to be very small. Our results can be summarized as follows: 

1. The differential cross section for solar neutrinos decreases sharply as the nuclear recoil energy 
increases. It almost vanishes beyond 1 keV (5 keV) for intermediate (light target), like ^^^Xe 
i'^^S). On the other the corresponding event rates for WIMPs of mass w 100 GeV extend 
further than 30 (150) keV ^^i^e (^25) respectively. 

2. The event rates for boron solar neutrinos at zero threshold energy and no quenching are 2- 
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FIG. 16: The ratio of the total cross section with threshold divided by that with zero threshold for A=131 
as a function of the threshold energy. Otherwise the notation is similar to that of Fig. 1121 Note that the 
observed events are an order of magnitude down, if the energy threshold is 1 keV. 
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FIG. 17: The ratio of the total cross section with threshold divided by that with zero threshold for A=32 
as a function of the threshold energy. Otherwise the notation is similar to that of Fig. 1161 Note that the 
observed events are an order of magnitude down, if the energy threshold is above 2 keV. 



3 orders of magnitude smaller than those for WIMPs with a nucleon cross section 10~^pb. 
Thus solar neutrinos are not a serious background down to 10^ ^"pb, but they may have to be 
considered at the level of 10^ ^^pb. 

3. Since the nuclear recoil energy arising from solar neutrino scattering is smaller than that asso- 
ciated with heavy WIMPs, one can further substantially decrease its contribution by restricting 
the observation of the recoil energy spectra above a few keV without seriously affecting the 
corresponding WIMP rates. Thus neutrinos do not appear to be a serious background even at 
the level of 10-"pb. 

4. By exploiting the quenching factors one may reduce this background still further. 



5. It should be noted that the solar neutrinos do not affect the DAMA result, Bernabei et al [2^, 
in both the energy and the cross section. DAMA uses the Nal target with a large quenching 
factor. Since events from Nal are mostly due to ^^^I, the event rate and the quenching factor 
are nearly the same as those for -'^^^Xe discussed in the text. Thus the solar neutrino events 
are below 1 keV, i.e. they are below the DAMA energy bin 2-4 keV. In other words the solar 
neutrinos may be dangerous for WIMP detection for WIMP-nuclcon cross section less than 
10~^° pb, which is far below the DAMA region of 10"^ - 10"'^ pb. 

6. The observation of the annual modulation of the signal (see e.g. [s^ and references there in) 
or even better by performing directional experiments [3^.[4o|. i.e experiments in which the 
direction of recoil is also measured, one will be able to select WIMP signals and discriminate 
against neutrino scattering. 

In the above discussion we focused on the coherent WIMP-nucleus scattering. We should not, of 
course, forget the spin contribution due to the axial current. In this case one has to deal with the 
proton and neutron spin nuclear matrix elements and the relevant elementary proton and neutron 
spin amplitudes. So the obtained results will depend on the specific target. It is clear, however, 
that the spin matrix elements do not exhibit coherence, i.e. do not scale with A^. Thus the event 
rates will be suppressed. In other words the boron neutrinos maybe be a serious background for 
nucleon spin cross section at the level of lO^'^pb. This, of course, will be worrisome, but then the 
corresponding coherent cross section must be less than lO^^^pb, since it is only then that the two 
modes can compete. 
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